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Abstract 

We consider all genus zero and genus one correlation functions for the Vi- 
rasoro vacuum descendants of a vertex operator algebra. These are described 
in terms of explicit generating functions that can be combinatorially expressed 
in terms of graph theory related to derangements in the genus zero case and to 
partial permutations in the genus one case. 
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1 Introduction 



A Vertex Operator Algebra (VOA) (e.g. [B], |FLMj . [K] ) is an algebraic structure 
closely related to Conformal Field Theory (CFT) in physics e.g. [BPZ] , [DMSJ. An 
essential ingredient of a VOA or CFT is the existence of a Virasoro vector whose vertex 
operator modes generate a Virasoro subalgebra of central charge C. The purpose of 
this paper is to describe all genus zero and genus one correlation functions for Virasoro 
descendants of the vacuum vector in terms of explicit generating functions. Examples 
of such Virasoro correlation functions have long been known in the VOA and CFT 
literature. For instance, they occur in describing the Kac determinant, CFT Ward 
identities and VOA modular differential equations e.g. |DMS] . [Z] . |MT2j . Our aim 
here is to provide a complete description of these Virasoro correlation functions at 
genus zero and one. 

We begin in Section |2] with a brief review of the theory of Vertex Operator Algebras 
(VOA). We describe a set of symmetric rational generating functions for all genus 
zero correlation functions for Virasoro vacuum descendants. These satisfy a recursion 
relation due to Zhu [Z] which in Theorem 12.31 imply that each generating function can 
be expressed as the sum of certain rational weights of appropriate graphs related to 
derangements in combinatorics. We also show that the generating function can be 
combinatorially expressed in terms of a /3-extension of a permanent |FZ] . 

In Section [3] we consider genus one correlation functions for Virasoro vacuum de- 
scendants. We begin with a review of elliptic and modular functions and prove in 
Theorem 13.31 that the Weierstrass function satisfies an important non-linear differen- 
tial equation. (As a by-product of this analysis we also obtain a closed formula for 
the modular derivative of any Eisenstein series). We then describe a set of symmet- 
ric elliptic generating functions for all such genus one Virasoro correlation functions 
which satisfy a genus one Zhu recursion relation [Zj. The resulting expression for the 
generating function is not manifestly symmetric. However, by use of the non-linear 
differential equation of Theorem 13.31 we show in Theorem 13.61 how to express each 
generating function as a sum of weights of appropriate graphs related to partial per- 
mutations. We also show that this generating function can be expressed in terms of an 
a, /3-extension of a partial permanent defined in Section I3.4[ a natural generalization 
to partial permutations of the notion of the /3-extension of a permanent defined for 
permutations. The Section concludes with an alternative graphical approach based 
on permutations. Finally, the Appendix describes the computation of the exponen- 
tial generating function for enumerating the various graphs employed in the earlier 
sections. 

2 Genus Zero Virasoro Correlation Functions 

2.1 Vertex operator algebras 

We review some basic properties of Vertex Operator Algebras e.g. [FLM] . |FHL] . 
H, [EE], [MS], [MT2] . A Vertex Operator Algebra (VOA) is a quadruple (V,Y,l,u) 
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consisting of a Z-graded complex vector space V = © ngZ V n with dim V n < oo, a linear 
map Y (Endy)[[z, z^ 1 ], for formal parameter z and a pair of states: the vacuum 
1 G Vo and a conformal vector u; G V2. For each u G V we have a vertex operator 

which satisfies the following axioms: 

Lower Truncation. For each pair u, v G V then v n u = for n ^> 0. 
Locality. (z x - ^ 2 ) iV [^(«, ^(u, z 2 )] = for N > 0. 
Creativity. Y(v, z)l = v + O(z). 

Virasoro Structure. There exists a conformal vector uj G V2 with 

y(w,z) = ^L n ^- n - 2 , 

where the modes L n = u n +i satisfy a Virasoro Algebra of central charge C with 
bracket relation 

n , _ , , (m 3 — rn) 

[■L> m i L n \ — [m — n)L m+n H — Oo m _ n . 

Furthermore, the Z-grading on V is provided by L with L v = nv for all v G V n . 
Translation. V(t; , z)] = 2). 

These axioms imply Y(l,z) = Idy and associativity: 

{z + z 2 ) M Y{u,z + z 2 )Y{v,z 2 )w = {z + z 2 ) M Y{Y{u,z )v,z 2 )w, (2.1) 

for u,v,w G V and M 3> [FHLj . In this paper we consider VOAs of CFT-type for 
which Vo = CI (and consequently, V n = for n < |MT2j ). 

2.2 Virasoro correlation functions 

Let v 1, • • • , v n G V and define the genus zero n-point correlation function for formal 
parameters Z\ , . . . , z n by 

(vi, zi; . . . ; v n , z n ) = (1', Y(vx, zi) . . . Y(v n , z n )l), (2.2) 

where V denotes the dual of the vacuum vector. Lower truncation and locality imply 
that the n-point function is a rational function in Zi, . . . , z n (expanded in the formal 
domain \zi\ > ... > \z n \). Thus it may be considered as a rational function on the 
complex Riemann sphere for Zi G C [FHL] . [MT2j . In particular, we consider the 
genus zero Virasoro m-point function 

G$(zi, ...,z m ) = Z { y ] (u, zi',,...;u, z m ) , (2.3) 

with Vi = . . . = v m = u, the Virasoro vector. We then find 
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Proposition 2.1. Gm(zx, • • • , z m ) is a symmetric rational function in Zi, . . . , z m . The 

set of such functions generate all genus zero n-point correlation functions for Virasoro 
vacuum descendants. 

Remark We refer to Gm(zi, ■ ■ ■ , z m ) as the order m genus zero Virasoro generating 
function. 

Proof Locality implies that Gm(zi, . . . , z m ) is a symmetric rational function in Z\, . . . , z n 
We illustrate the generating property by considering the 2-point function for two Vi- 
rasoro vacuum descendants V\ = L_i x . . . L_i r l and v 2 = L-fo . . . L_/. s l. Consider 

(1', Y(Y(u, wt)... Y(co, w r )l, w)Y(Y(u, z x ) . . . Y(u, z.)l, z)l). (2.4) 

The two point function for v\, V2 is the coefficient of Yli=i YYj=i z / 1 i* 1 d23|). By 
repeatedly using associativity ( 12. D we find 



l[(w r + w) M Y[{ Zj + z) N (l', Y(Y{u>, Wl )... Y{lo, w r )l, w)Y(Y(lo, Zl ) . . . Y(u, z s )l, z)l) 

i=l j=l 
r s 

= Y[(wi + w) M Y[( Zj + z) N (l', Y(co, wi + w)Y(Y(oo, w 2 ) . . . Y(co, w m )l, w) 

i=l j=l 

Y(cu, Zl + z)Y(y(w, z 2 )... Y(cu, z n )l, z)l) 

r s 

= (1', Y[( Wl + w) M Y[( Zj + z) N Y(u, Wi + w)Y(u, Zj + z)l) 

i=l j=l 
r s 

= Y\.( Wi + II + Z ) N ^r+s( w l +W, . . . ,W r + W,Zx + Z, . . . ,Z S + z), 

i=l j=l 

for M, N ^> 0. Thus ( I2.4p and g!£+ s (wi+w , . . . , w r +w, Zi+z, . . . , z s +z) are determined 

by the same rational function whose coefficient of YYi=i ^i" 1 Y\]=i z j J 1 i s the 2-point 
function for Vi, v 2 . The general result follows from a similar argument. □ 

By Lemma 2.2.1 of [Z] we obtain a recursion formula for the order n Virasoro 
generating function: 

Proposition 2.2. 

'AZ 2 , ■ ■ . , z r> 



V lfe 



C 1 

+ -77—G n _ 2 (z 2 , . . . , z k , . . . , z n ), (2.5) 
A z ik 

where Zij = Z{ — Zj and z^ denotes that z k is omitted. □ 

Thus we find using G^ = 1 and Gf'\zi) = that 

n( Q )< \ CI ^(o)/ ^ C 

(^2 l z l> Z V — 7T ^3 \ z l, z 2,Zz 



9 _4 ' "3 V~±>~^~.V 2 r 2 r 2 • 

z ^12 ^12^13 z 23 
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2.3 Genus zero Virasoro graphs 

We next develop a graphical/combinatorial approach for computing the order n Vi- 
rasoro generating function. To this end we define an Order n Genus Zero Virasoro 
Graph to be a directed graph Q with n vertices labelled by z\ , . . . , z n with the following 
properties. Each vertex has degree 2 with unit indegree and outdegree and no vertex 
is joined to itself. Thus the connected subgraphs of Q are r-cycles with r > 2 vertices. 
We define a weight W(Q) on Q as follows. For each directed edge with i ^ j, 

Zi O K3 Zj 

define an edge weight 

W{£, 13 ) = -L. 

Let K be the number of cycles in Q and define 

K 



where the product ranges over all the edges of Q. 

In the Appendix we describe a generating function for counting the number d Ut K 
of inequivalent order n graphs with K cycles. This is determined by the polynomial 

dM = £ d n , K p K = (-l)»n! £ ^ ( n "_f <) , (2-6) 
x>i i=o ' ^ ' 

for cycle counting parameter /3. Thus for n = 2 we find d 2 (P) = (3 so that there is one 
Virasoro graph Q with weight 

the order 2 Virasoro generating function. For n = 3 we have d 3 (/3) = 2/? so that there 
are two order 3 Virasoro graphs: 



zi 





Gl ~ Z 3 d^-t)Z 2 > & - ^ ^ 

The sum of the weights of these two graphs is 

w(g 1 ) + w(g 2 ) = f -r4-r + ?-r^r = gP(z 1 ,z 2 , z 3 ), 

z ^12^23^31 z ^13^32^21 

the order 3 Virasoro generating function. For n = 4 one finds d A (f3) = 3/3 2 + 6/3 so 
that there are 9 independent graph configurations whose weights sum to the order 4 
Virasoro generating function. These examples illustrate the following general result: 
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Theorem 2.3. The order n genus zero Virasoro generating function Gn^ is the sum 
of the weights of all inequivalent order n genus zero Virasoro graphs. 

Proof We prove the result by induction in n. The result holds for n = 2 and 3 as 
described above. For n > 4 we reinterpret the recursion formula equation (12. 5p in 
terms of a construction of order n Virasoro graphs from order n — 1 and n — 2 Virasoro 
graphs. 

Divide the set of inequivalent order n Virasoro graphs into the following two types: 

Type I. Graphs for which the Z\ vertex belongs to a 2-cycle. 

Type II. Graphs for which the Z\ vertex belongs to an r-cycle for r > 3. 

Consider the following term arising in the RHS of (12 .5p for k — 2, . . . , n 

CI 
~2~z 



4 Gn-2( z 2> ■ ■ ■ > *k, ■ ■ ■ , Z n ). (2-7) 



life 

By induction, G^l 2 i s ^ ne sum °f weights of inequivalent Virasoro order n — 2 graphs 
j w ith vertices labelled by Z2, ■ ■ ■ , ifc, ■ ■ ■ , z n . For each such Q n ~ 2 , let denote 
the order n Virasoro graph formed by the product of Q n ~ 2 and the 2-cycle with Z\ and 
z k -vertices. Clearly, is an order n graph of Type I and all inequivalent such graphs 
arise in this way. Since 



C 1 



then summing (12. 7p over all k results in 



E ^ G n-2(Z2, ...,Z k ,...,Z n ) = J2 W(W), {2.1 



k>2 Jlk g™ 
where ranges over the inequivalent Virasoro order n graphs of Type I. 

We next consider the remaining terms arising on the RHS of (12. 5 j) for k — 2, . . . , n: 

T^G ( ^}_ l (z2, ■ ■ ■ , z n ), 

where V = ^^=2(7"^ + ~£~)- By induction, G^l 1 (z 2 , ■ ■ ■ , z n ) is the sum of weights 

1^ Ife 

of all inequivalent order n — 1 Virasoro graphs Q n ~ l labelled by Z2, ... , z n . Let {Sij} 
denote the edges for such an order n — 1 graph. Recalling that W(£ij) = \ we find 

that 

VW{G^) = W(G^)zl (^d Zi + ±d Zj + -L + -L) * 

1~ \zu ztj z u z Xi ) z ii 



Using the differential equation 



-d x --d v + — + -\ = — , 2.9 

x y x z y z J (x — y) z x z y z 



for x = Zu and y = z\j we obtain 



where Gfj denotes the order n Virasoro graph obtained by "inserting" a vertex labelled 
Z\ between each connected pair and Zj in Q n ~ x i.e. schematically 

Zi 

■ ■ • Zi O KD Zj ... >. . . . Zi O >0 K3 % 

Clearly, each order n Virasoro graph Qf- is of Type II. Summing over all inequivalent 
order n — 1 graphs {{? n_1 } results in a sum over all inequivalent graphs {Q\i\ of Type II 
i.e. 

y ll 

Combined with the Type I graphs from (12.81) we therefore find 

Gfte, . . . , z n ) = Y, W(G?) + Y W(0n), 

as required. □ 



2.4 Derangements and the /^-extension of a permanent 

The set of order n Virasoro graphs {Q} is in 1-1 correspondence with the set of de- 
rangements (i.e. fixed point free permutations) of the label set $ = {1, . . . , n}. 
Thus the directed edges {£ij} of Q correspond to a derangement permutation cf> e F(<&) 
with <fi(i) = j i. In particular, connected cycle subgraphs correspond to cycles of 
the corresponding permutation. Hence we may write the result of Theorem 12.31 as 

/ r \K{<t>) n , 
<j>EF($) ^ ' i=l i( t>(ki) 

where K(<p) is the number of cycles in the derangement permutation 0. The expres- 
sion on the RHS of (12.101) can be expressed in terms the so-called /3 -extension of a 
permanent of an n x n matrix A = (a^) defined by [FZj 

n 

pernig A = £ Y[a ln{i) , (2.11) 

7TgS n i=l 

where 7r ranges over all permutations S n of the set {1, . . . , n) and K(tt) is the number 
of cycles in n. (For j3 = 1 then perm x A = perm A, the standard permanent, whereas 
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for (3 = — 1 we obtain perm_ 1 A = (— l) n det^4). In the present context, we define a 
symmetric matrix A with components 



(hj = s i 



i = j, 



Since the diagonal elements of A are 0, only fixed point free permutations contribute 
to (12. lip so that Theorem 12.31 can be restated as 

Theorem 2.4. The order n genus zero Virasoro generating function is given by 

G^\zi, . . . , z n ) = permc A. □ 



3 Genus One Virasoro Correlation Functions 
3.1 Elliptic functions and modular forms 

We first consider aspects of the theory of elliptic and modular functions relevant to 
our discussion. We note the following conventions: q z = e z , q = q 2 ^ T = e 2mT and 
d = qd q = i^pdr for z e C and r € H, the complex upper half plane. Define 

P 2 {z,q) = p( Zl T) + E 2 (q) 

1 oo 

= -a+X)^- 1 )^?)^" 2 . (3- 1 ) 

k=2 

where p(z,r) is the Weierstrass function periodic in lui and 2-iriT e.g. [WW] . [C] . 
Ek(q) is zero for k odd, and for k even is the Eisenstein series 

v ' n>l y 

where B^ is the Bernoulli number determined by 

1 ' k\ 



q z - 1 ^ k\ 2 12 



Ek(q) is a holomorphic modular form of weight k on SX(2,Z) for > 4 [S] whereas 
E 2 (q) is a quasi- modular form [KZ] . We also define Pk(z, q) for k > by 



P (*, 9) = - log(«) + X \ E M*\ (3-3) 
9) = -j^fciW*. q) = ~ k + (-1)' X ( ]fc _ ! ^ ( 3 - 4 ) 
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Pi is periodic in 2iri but is not periodic in 1-kit but P k is periodic in both 2ni and 
2-kit for all k > 2. Pq(z) is related to the elliptic prime form K(z, q) |Muj 

K(z, q) = ex P (-P (z, q)) = (3-5) 

where 



n>l 



are the Jacobi theta series and the Dedekind ^-function respectively. 
We find the following: 

Proposition 3.1. 

2dP (z, q) = P 2 (z, q) - P?(z, q) - 3E 2 (q). (3.8) 
Proof We first note from ( 13T2j) and (J3~Z3) that 

dlog( V (q)) = ~E 2 (q). (3.9) 

Furthermore, 9i(z,q) obeys the heat equation 

2d6 1 (z,q) = d 2 z 9 1 (z,q). (3.10) 

It thus follows from (13.5P that 

2dd x (z,q) = (-2dP (z 1 q)-3E 2 (q))9 1 (z,q), 
d%(z,q) = (Pl(z,q)-P 2 (z,q))e 1 (z,q). 

Hence the heat equation (13.101) implies the result. □ 

Comparing the Laurent expansions in z of both sides of ( 13. 8 p immediately implies 

Proposition 3.2. The Eisenstein series obeys 



k(k + 3) k 



k 



dE k {q) = ^ r ^E k+2 {q)--Y J EME k+2 _ r {q)- □ (3.11) 



r>2 



Remark Let D k = d + kE 2 (q) denote the modular derivative which when applied to 
a modular form /(r) of weight k results in a modular form D k f(r) of weight k 4- 2 [Z]. 
|Maj . For k > 4, Proposition 13.21 implies 



D fc £ fe (g) = k ( k + 3) E k+2 (q) - ~ J] E r (q)E k+2 „ r (q), 



r>4 



i.e. the weight + 2 modular form D k E k (q) is explicitly known in this case. 
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We now consider a non-linear differential equation satisfied by the Weierstrass 
function analogous to (12. 9p for rational functions. 

Theorem 3.3. 

(d- Pi(x,q)d x - P 1 (y,q)dy + P 2 (x,q) + P 2 (y,qf)P 2 (x -y,q) = P 2 (x , q) P 2 (y , q) . 

(3.12) 

Proof Apply d\ to fl3j| to find dP 2 (z,q) = 3P A (z,q) - P 2 (z,q) - 2P 1 (z, q)P 3 (z, q). 
We then find that the LHS of (13.121) is symmetric and elliptic in x, y with poles only 
at x = and y = 0. The result follows on comparing the leading terms of the Laurant 
expansion in the neighbourhood of the poles on both sides of (13.121) . □ 

We also note the identity 

dE 2 (q) + 2E 2 (q)P 2 (x, q) + P 4 (x, q) = P*(x, q), (3.13) 

obtained from (13.12p in the x — > y limit. 

3.2 Genus one n-point correlation functions 

In order to define and compute n-point correlation functions at genus one, Zhu [Z] 
introduced the so-called square bracket formalism. This consists of a second VOA 
(V, Y[ , ], l,d)) isomorphic to (V,Y( , ), l,oo) with vertex operators 

Y[v,z] = ^M^' 1 = Y{q L z °v,q z - 1). 
The new conformal vector is to = u — ^1 with vertex operator 

The genus one partition (or 0-point) function for V is defined by 

Z^\q) = Tr v (q L °- c / 2i ) = ^ dim V n q n ~ c ' 24 . (3.14) 

n>0 

For v\ , . . . , v n 6 V we define the genus one n-point correlation function for parameters 
z 1 ,...,z n by 

4 X) (vi, zi; . . . ; v n , z n , q) = Tr v {Y(q%>v u q 9l ) . . . Y(q^v n , q Zn )q L °- C/M ) . (3.15) 

Zhu describes a recursion formula for expanding an n-point correlation function as a 
linear combinations of n— 1 point correlation functions with coefficients given by elliptic 
functions [Zj. Thus genus one n-point correlation functions are elliptic functions just 
as genus zero n-point correlation functions are rational functions. 

We next define the order m genus one Virasoro generating function by 
G^{zx, ...,z m ;q) = Z { y ] (u, zi;...;u, z m ; q) . 
Analogously to Proposition 12.11 we then find, using Lemma 5 of |MTlj , that 
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Proposition 3.4. (z b . . . , z m ; q) is a symmetric elliptic function in Zx,...,z m . 
The set of such functions generates all genus one n-point correlation functions for 
Virasoro vacuum descendants. 

Applying the Zhu recursion formula (Proposition 4.3.2, [Z]) we obtain the following 



genus one analogue of Proposition 12.21 

Proposition 3.5. The genus one Virasoro generating function obeys the recursion 
relation 

GS\z!, z 2 , ... , z n ; q) = dG^fa, . . . , z n ; q) 

n 

+ Pl ( Zlk > <l) d Zk G n-l( Z 2, ■■■,Zn,q) 



k=2 

£>ri\&Xki q)^ 

k=2 



q)G^ ) _ 1 (z 2 ,...,z n ;q) 



+ Y ^"^(^lfc, q)G^l 2 (z 2 , . . . , z k , . . . , z n ; q), (3.16) 

k=2 

where £k denotes that z k is omitted. □ 

The RHS of (I3.16P is not manifestly symmetric in z±, . . . , z n as required by Propo- 
sition 13.41 We show below in Theorem 13.61 how to express Gn in a symmetric fashion 
as a sum of weights of appropriate graphs. In order to achieve this it is useful to define 

©v(?) = v c (q)4\qI 

T n (z 1 ,...,z n ;q) = rf ' {q)G l ^ ) (z u z 2 , . . . , z n ; q), 
for Dedekind rj(q) of (13. 7p . Thus T (q) = 0y(<?)- Applying (I3.9P we find 

T 1 (z 1 ;q)= (d+jEtiqfj Q v {q). 

For n = 2 we apply (13.131) to obtain 

I 2( '-I- | 0~ ' CE,( U )0- [ ^-E 2 (q)) + 2P 2 (z 12 , q)d + %P 2 2 {z X2 , q) 1 Q v (q). 
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For n = 3 we apply (13.121) to eventually obtain 




+ —E 2 {q) (P*(z 12 , q) + P*(z 13 , q) + P 2 2 (z 23 , q)) 
+ 2 (P 2 (z 12 , q)P 2 (z 13 , q) + P 2 (z 12 , q)P 2 {z 23 , q) + P 2 (z 13 , q)P 2 (z 23 , q)) d 



+ CP 2 {z 12 , q)P 2 {z 13 , q)P 2 (z 23 , q) Q v (q). 



(3.17) 



This expression depends on C, E 2 (q), P 2 (zij, q) and Qv{q) and is manifestly symmetric 
in z 1 ,z 2 ,z 3 as expected from Proposition 13.41 

3.3 Genus one Virasoro graphs 

Similarly to Section I2T31 we now develop a graphical/combinatorial approach for com- 
puting the order n genus one Virasoro generating function. We define an Order n 
Genus One Virasoro Graph to be directed a graph Q with n vertices labelled by 
Z\, . . . , z n . Each ^-vertex has degree deg(zj) = 0,1 or 2. The degree 1 vertices can 
have either unit indegree or outdegree whereas the degree 2 vertices have both unit 
indegree and outdegree. Thus, in this case, the connected subgraphs of Q consist of 
r-cycles for r > 1 vertices and necklaces with degree 1 end-vertices and possible inter- 
mediate degree 2 vertices. We regard a single degree vertex as being a degenerate 
necklace. 

We define a weight W(Q) on Q as follows. For each directed edge £y define an edge 
weight 



Let K be the number of cycle graphs in Q and let M be the number of necklaces and 
define 



where the product ranges over all the edges of Q. Thus the weight depends on C, 
E 2 {q), P 2 (zij,q) and Q v (q)- 

In the Appendix we describe a generating function for counting the number p n ,K,M 
of inequivalent order n Virasoro graphs containing K cycles and M necklaces deter- 




(3.18) 




(3.19) 
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mined by the polynomial 

Pn(a,P)= T Pu,k,m^ M P K = {~l) n n\ V ^ ~ A (3.20) 
' z! \ n — i J 

K>0,M>0 i=0 v 7 



for necklace and cycle counting parameters a and j3 respectively. Thus for n = 1 we 
find pi(a, (3) = a + (3 corresponding to the two inequivalent graphs 

Gi = z io, Q 2 = Cp z ± , 

i.e. a degenerate necklace and a 1-cycle. Clearly the sum of their weights is 

W(Gx) + W(G 2 ) = d& v (q) + jE 2 (q)Q v (q) = T^z,; q). 

For n = 2 we have p 2 (a, j3) = a 2 + 2a(3 + (3 2 + 2a + (3 for the 7 inequivalent graphs: 



'i a O z 2 , Zi O 



ZiO O ^2 ' Zl " ° Z<1 ! Zl ° O 22 J 21 O > 

2 lC) ^ 0^2' Zl ° *°*2, Z lO Z2 - 

whose weights sum to 

YW{Q) = d 2 Q v (q) + ^E 2 (q)dQ v (q) + ^E 2 (q)dQ v (q) + ^-E 2 (q) 2 Q v (q) 

g 

+ P 2 (z 12 , q)dQ v (q) + P 2 (z 21 , q)de v (q) + ^P%(z 12 , q)G v (q) 
= T 2 (z 1 ,z 2 ;q). 

For n = 3 one finds 

p 3 (a, f3) = a 3 + 3a 2 f3 + 6a 2 + 3a/? 2 + 9a/3 + 6a + /3 3 + 3/3 2 + 2/?, 

describing 34 inequivalent genus one Virasoro graphs whose weights sum to Ts(zi, z 2 , z^; q). 
These examples illustrate the general result: 

Theorem 3.6. T n (zi, . . . , z n ; q) is the sum of weights of all inequivalent order n genus 
one Virasoro graphs. 

Proof We prove the result by induction in n. We have already shown the result holds 
for n = 1 and n = 2. For n > 3 we reinterpret the Zhu recursion formula (13.161) in 
terms of a construction of order n genus one Virasoro graphs in terms of order n — 1 
and n — 2 graphs. 

Let Q n denote an order n Virasoro graph. Based on the nature of the zi-vertex, 
this must belong to one of the following five types: 
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Type I. degOO = 0: 

ZiO 

Type II. deg(zi) = 1: 

Z\ o K) z a ■ ■ ■ or zi O o z a - • • 

Type III. deg(zi) = 2 where the zi-vertex forms a 1-cycle: 

-•.O ■•• 

Type IV. deg(^i) = 2 where the zi-vertex is an element of a 2-cycle: 

ZlO^Z^Za ■ ■ ■ 

Type V. deg(zi) = 2 where the zi-vertex is an element of a necklace or an r-cycle 
for r > 3: 

Z\ 

• ■ ■ Z a O KD >0 Zb ■■■ 

Applying the Zhu reduction formula (13 . 16[) we have 

V n {z u ...,z n ) = rfd [vr c Y n ^{z 2 , z n )) (3.21) 

n 

+ p i( z ik)d Zk Tn-i{z2, ...,z n ) (3.22) 



k=2 
n 



+ J2^2(zik)T n - 1 (z 2 ,...,z n ) (3.23) 



k=2 
n 



+ ^2 -^ p i( z ik)^n-2(z 2 , . . . , z k , . . . , z n ). (3.24) 

k=2 

Here, and below, we suppress the (/-dependence as a convenient abbreviation. 

We now show how the different parts of Zhu's reduction formula relate to the 
Virasoro graph weights by using induction to express both T^ and r^_ 2 as the sum 
of weights of all order n — 1 and n — 2 Virasoro graphs, respectively. Let Q n ~ l denote 
an order n — 1 Virasoro graphs labelled by z 2 , ■ ■ ■ , z n . Then by induction we have 

r n _! = r n ^(z 2 , ...,z n ) = J2 wig*- 1 ), 

gn-l 

for 

W{Q n ~ l ) = (^) K l[w(S)d M Q v , 

£ 

where Q n ~ x consists of M necklaces and K cycles with edges {£}. 
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Consider the first contribution (I3.2ip arising from Zhu reduction. Recalling (13. 9p 

c 

2 " 

K 



then d(r]- c Q v ) = r]~ c (^E 2 + d)& v so that 



rfd ( V ~ C W (G^)) = (|) m W (£)) d M ^Q v (3.25) 

+ E 2 ^ ^ W (S)\ d M e v (3.26) 

+ (f) i "f a n w (Oj^©v- (3-27) 

The first term (I3.25P is the weight of the order n Virasoro graph £/" of Type I consisting 
of a Zi-labelled degenerate necklace together with Q n ~ l . Clearly, all inequivalent Type I 
graphs can arise in this way. The second term (I3.26P is the weight of the order n 
Virasoro graph Q^j of Type III consisting of a Zi-labelled 1-cycle together with Q n ~ l . 
Likewise, all inequivalent Type III graphs arise in this way. Thus summing (I3.25P and 
(EEggD over all Q n ~ x we find 

E (Tf 11^) (° M+1 + ^ E * dM ) Qv = J2 W ^ + J2 W ^), (3.28) 

summing over all inequivalent order n graphs of Type I and III. The remaining term 
(I3.27P contributes 

(§) {d\{w{£)^d M e v = w{g n - l )Y^d\o g w{£). (3.29) 

In a similar fashion, (I3.22|) of the Zhu reduction formula gives a Q n ~ l contribution of 

n 

W{G n - v ) E P ^)d Zk log W(£). (3.30) 

e k=2 

Likewise, (13. 23j) gives a Q n ~ x contribution of 

w{G n ~ x ) J2 m{zi k ) = w{G n - x ) I ( p ^) + P2{z lb )) + J2 P2 (^) J > 

k=2 \£ab c,deg(z c )=l J 

(3-31) 

where the first sum is taken over all Q n 1 edges £ a b and the second sum over all 
^-vertices of degree 1 arising from necklace end-vertices. Let 



£ah 



(d + J2 p i( z ik)d Zk ) logI^(£ afe ) + P 2 (z la ) + P 2 (z lb ) 

k=2 

+ p 2(^i c )), (3.32) 

c, deg(z c )=l 
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denote the sum of (I3.29l) - (l3.3ip . Consider an edge £ ab with z a , z&-vertices contributing 
to (I3.32p where either (A) z a = z b or (B) z a ^ z b . 

Case (A). Here the z a -vertex forms a 1-cycle with edge weight W(£ aa ) = E 2 . 
Hence W(Q n ~ 1 ) = ^E 2 W(Q n ~ 2 ) where Q n ~ 2 is an order n — 2 graph labelled by 
z 2 , ■ ■ ■ , z a , . . . , z n (with label z a deleted). All distinct such order n — 2 graphs arise in 
this way. The £ aa dependent part of (I3.32p is thus 

^W(G n ~ 2 )E 2 (d\ogE 2 + 2P 2 (z la )) . (3.33) 

The last term (I3.24p in the Zhu reduction formula makes a Q n ~ 2 contribution of 
<2 W(G n ~ 2 )P 4 (z la ). Adding this to (13331) and using f[3TTB"j) we obtain 

W(G n - 2 )j(dE 2 + 2E 2 P 2 (z la ) + P 4 (z la )) = W(G n ~ 2 )jP 2 (z la ) 2 



IV J J 



where C/fy is the order n graph of Type IV formed by the union of Q n 2 with the 
2-cycle with vertices labelled by z\ and z a . Summing over all z a and Q n ~ 2 we obtain 

Case (B). In this case W(S a b) = P2(z a b) so that the total E ab contribution arising 
from the edge sum in S{Q n ~ l ) of (13.321) is 

W(^ n - 1 )P 2 (2 o6 )- 1 (d + P^d^ + P 1 (z lb )d Zb + P 2 (z la ) + P 2 {z lb ))P 2 {z ab ) 

= W{0 n-l ^a)Hz lb ) ^ (3 
P2\Z a b) 

using Theorem 13.31 with x = —z la ,y = —z u . But the RHS of (13.341) is the weight of 
the order n graph Qy of type V constructed by inserting z\ between the z a -vertex and 
the z b -vertex: 

■■■ Z a O KD Z b ... 

Summing over inequivalent Q n ~ l we obtain Y2g n W^(^v), summing over all inequivalent 
graphs of type V. 

Lastly, if deg z a — 1 then the z a -vertex is a necklace end-vertex and we obtain an 
additional contribution from the last term in (I3.32p of 

w(g n - l )p 2 (z la ) = w(g? 1 ), 

where C?fj is the order n graph of Type II found by joining the zi-vertex to the z a - 
vertex. Summing over inequivalent Q n ~ x and all such necklace end z a -vertices results 
in ^2gn W(Qij). Thus altogether we have shown that 

T n (z 1 ,... 1 z n ,q) = J2W(g n ), 

gn 

summing over all inequivalent order n Virasoro graphs of type I-V and hence the 
theorem holds. □ 
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Remark Theorem 13.61 can be readily generalized for any ^/-module M by replacing 
®v(q) by Qm(q) — r)(q) c Tr M (g^ ) -0 / 24 ) throughout. The remaining dependence on 
C, E 2 (q) and P 2 (zij,q) remains unchanged. 



3.4 Partial permutations 

Let ty n denote the set of partial permutations of the label set $ = {1, . . . , n} i.e. the 
injective partial mappings (including the empty map) from {1, . . . ,n} to itself. The 
set of order n genus one Virasoro graphs is in 1-1 correspondence with \I/ n where the 
directed edges of Q correspond to ip G ^ n with = j. 

It is natural to introduce the Partial Permanent of an n x n matrix A = (ay) 
indexed by i, j G {1, . . . , n} as follows 

pperm A = ^ | J a^ w , (3.35) 

i/)gr J6 dom ^ 

where dom ijj denotes the domain of the partial map ifi and with unit contribution for 
the empty map. We also define the a, /3 -Extended Partial Permanent by 

Wenn af) A=Y,<* NW P KW U (3-36) 

where N(ip) is the number of necklaces and K(ip) is the number of cycles determined 
by i>. 

We now utilize partial permutations to re-express the order n genus one Virasoro 
generating function. Towards this end, we define a linear map C a from C[a], the vector 
space of complex coefficient polynomials in a, to the complex vector space spanned by 
d k ®v(q) as follows 

n n 

C a (Y^c k a h ) = Y,c k d k Q v (q), c k G C. (3.37) 

k=0 k=0 

Then Theorem 13.61 can be restated as 

Theorem 3.7. The order n genus one Virasoro generating function is given by 

T n (zx, ...,z n ;q)=C a (pperm Q c , □ 
with W = (W lj ) for edge weight W {j of (ETISjl . 
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3.5 Further combinatorics 



Lastly, we briefly consider an alternative way of computing T n (zi, . . . ,z n ;q) in terms 
of a sum over the permutations of the label set $ = {1, . . . , n} by means of an alter- 
native graphical approach. Define a truncation operation T p on the vector space of 
polynomials C[p] by 

n 

T p (22 c kP k ) = c o + cip, c k e C. 

k=0 

We also define a linear map C p as in (I3.37p . 

We next define an Order n Genus One Permutation Virasoro Graph to be a directed 
graph Q with n vertices labelled by Z\, . . . , z n where each vertex has degree 2 with unit 
indegree and outdegree. The set of such graphs is in one to one correspondence with 
the permutation group S n on n symbols. We define a weight x{G) on Q as follows. 
For each directed edge £ij define a weight 



Xij 



^p + E 2 (q) i = j, 
^p + P 2 (z ij ,q) i^j. 



Then for a cycle subgraph Q a with edges £ a we define: 

C 

x{Gv) = -wT p ([[xij)- 
Each a corresponds to permutation group cycle. Finally we define the weight of Q by 



For example, there are 3! order 3 permutation graphs 



zic/ dz 2 dz 3 , g 2 = z 1 a^yjz 2 cfz 3 , fe= ^cO 3 ^ oz 2 



Qi = z 2 o^£)Zz cfz 1} g, 



5 




Z36* Oz 2 ' ^° z 2 6* OZ 3 




in contrast to the 34 order 3 graphs arising in the previous section. The sum of the 
weights for the displayed 6 graphs results in T 3 (zi, z 2 , z 3 ; q) of (J37TTJ) . In particular, 
x(Qi) gives the first line, x(£? 2 ) + x{<3z) + x(Ga) gives l ines and x{<3s) + x(£e) gives 
lines 6-7 of (13.171) . This illustrates the general result: 

Theorem 3.8. Tj, (zx, . . . ,z n ;q) is the sum of weights of all inequivalent order n genus 
one permutation Virasoro graphs. □ 
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The proof proceeds along similar lines to Theorem 13.61 Using the correspondence 
with E n we can also restate this result as a sum over permutations: 

r n (z u ...,z n ;q) =£ p | (2" ) n^fn^w))' ( 3 ' 38 ) 
where tv has cycle decomposition % = o\ . . . ok with cycles = (...«...). 



A Counting Graphs 

A.l Counting genus zero Virasoro graphs 

Let d nj K denote the number of order n genus zero Virasoro graphs with K cycles as 
introduced in Section [231 Define a polynomial d n ((3) = J2k>i d n ,Kfi K with d ((3) = 1. 
Recall that each order n genus zero graph can be expressed as a derangement on 
n symbols. Hence d n (l) is the number of derangements on n symbols. Define an 
exponential generating function for genus zero Virasoro graphs 



n>0 

We then find 
Proposition A.l. 

D(J),z)=(j^ • (A.l) 

Proof Define a weight function wp >z on each order n graph with K cycles by 

w PtS {G) = P K z n . 
The exponential generating function is thus 

g n - 

where n denotes the order of Q. Each genus zero Virasoro graph Q is the labelled 
product of r-cycles for r > 2. But wp iZ is multiplicative with respect to this cycle 
decomposition with Wp jZ {Q a ) = j3z r for any r-cycle. Thus we find e.g. |FSj 

D{P,z) = exp(C{P,z)), 

where C(/3, z) = J2g a ^iVp, z {Go) is the exponential generating function for all r-cycles 
for r > 2. Since there are (r — 1)! inequivalent r-cycles this is given by 



r>2 

Therefore the result follows. □ 
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It is straightforward to expand flA.lj) to obtain ( 12. 6 p and the following recursion 
formula 

dn{P) = (n - l)(4-i(/3) + Pdn- 2 (P)), = 0, do(/3) = 1. 

These agree with standard results for derangements in the case (3 = 1 e.g. [FSJ. 

A. 2 Counting genus one Virasoro graphs 

Let p n ,K,M denote the number of order n genus one Virasoro graphs with K cy- 
cles and M necklaces as introduced in Section [3.31 Define a polynomial p n (ot, 0) = 
^2k>o m>o Pn,K,M aM with po(a, 0) = 1. Define the exponential generating function 
for genus one Virasoro graphs 



~7J 

P(a, P,z) = y2 —rPn(oc, 0). 

— ' n.i 



We then find 
Proposition A. 2. 



n>0 



P(a,P,z) = jf^. (A.2) 



Proof Here we define a weight function w a ^ >z on each order n graph with K cycles 
and M necklaces by 

WaMG) = a M p K z n , 

so that 

P(a, P,z) = Y] —rW aAz {G), 



Q 



summing over all genus one Virasoro graphs Q for all order n. But the weight function 
is multiplicative with respect to the necklace Q n and cycle Q a decomposition so that 
here we find e.g. |FS] 

P(a, p, z) = exp(N(a, z) + C(p, z)), 

where N(a, z) = J2g N ^ w o,i3,z(Gn) is the exponential generating function for connected 
necklaces and C(P,z) = J2g a ^\ w a,p,z{Ga) is the exponential generating function for 
connected r-cycles for r > 1. Thus, much as before 

C(P, z) = pJ2 \/^ r " !) ! = "0 Ml - *)• 

Similarly, w a ,p, z (GN) = oiz r for a necklace with r > 1 vertices and since there are r! 
inequivalent such necklaces we find 



z 

az' r! = Ct- 



rl 1 — z 

r>l 

Therefore the result follows. □ 
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flA.2j) implies (I3.20p and the following recursion formula 



p n +i(a, P) = (2n + a + f3)p n (a, (3) - n(n + (3 - l)p n _ 1 (a, (3), 
with p-i(a, /?) = and po(c>o /?) = 1- 
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